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Abstract 

The energy spectrum in three examples of inhomogeneous, anisotropic turbulence, namely, 
purely mechanical wall turbulence, the Bolgiano-Obukhov cascade and helical turbulence, is an- 
alyzed. As one could expect, simple dimensional reasoning leads to incorrect results and must 
be supplemented by informations on the dynamics. In the case of wall turbulence, an hypothesis 
of Kolmogorov cascade, starting locally from the gradients in the mean flow, produces an energy 
spectrum which obeys the standard k~3 law only for kx% > 1, with x 3 the distance from the 
wall, and an inverse power law for kx% < 1. An analysis of the energy budget for turbulence in 
stratified flows, shows the unrealizability of an asymptotic Bolgiano scaling. Simulation with a 
GOY model, leads instead to a k~ a spectrum for both temperature and velocity, with a ~ 2, 
and a cross-correlation between the two vanishing at large scales. In the case of not reflection 
invariant turbulence, closure analysis suggests that a purely helical cascade, associated with a k~"s 
energy spectrum cannot take place, unless external forcing terms are present at all scales in the 
Navier-Stokes equation. 

PACS numbers: 47.27.-i, 47.27.Nz, 47.27.Te 
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I. Introduction 



Turbulence in nature is always inhomogeneous; the reason is the origin of the fluctuations, 
either in the instability of flow patterns, or in the presence in some finite volume, of tempera- 
ture gradients, chemical reactions or external stirring. If the Reynolds number is large, however, 
there are turbulent fluctuations at scales much smaller than that of the forcing, and to them, the 
idealization known as homogeneous isotropic turbulence can be applied. 

In practical applications, what one is interested in, is the effect of turbulence on the mean 
flow and on transport, which is parametrized in terms of eddy viscosities and diffusivities (see 
e.g. [[!]). In some cases, in order to calculate these quantities, some information on the turbulent 
energy spectrum is necessary, and Kolmogorov scaling is usually assumed. For instance, in 
the derivation of Lagrangean diffusion models J|, || , the Kolmogorov scaling hypothesis is present 
explicitly through the assumption of Markovian velocity increments, at time-scales below that of 
the energy containing eddies. 

What becomes necessary then, is some matching condition at the transition from the inertial 
range (small scale) to the energy containing range (forcing scale) ; this is essentially the problem of 
connecting a region of k~a scaling, to the peak in the energy spectrum. Phcnomcnological theories 
have dealt with this problems H ; more recently, the question of how fast the effect of inhomogeneity 
decays at small scales has been put under exam both theoretically and using reduced models |(| (t| . 

There are situations, however, in which the problem of how far the inertial range preserves 
memory of the inhomogeneity of the forcing, becomes particularly serious. 

The most obvious way this can happen is when forcing takes place at all scales. Notice that 
this does not necessarily require the presence, to make an example, of obstacles of corresponding 
sizes in the flow. Already in the case of wall turbulence || , one has mean flow gradients at lengths 
ranging from the viscous range to the height of the boundary layer, which leads to a situation of 
coexisting "energy range" and "inertial range" eddies, distributed at all scales. 

An extended forcing range develops clearly, also in the presence of stratification, due to the 
effect of buoyancy. In this case an additional scale, the Obukhov length marking the transition 
from mainly mechanical to convection dominated turbulence, becomes important, and the way 
in which mechanical and convective contributions to the dynamics balance one another, makes a 
description based on scaling rather non-trivial. 

A third way, in which the small scale dynamics of turbulence could be modified by processes in 
the energy range, is when helicity is fed, together with energy, into the system. The importance 
of this process has been discussed recently by Yakhot (7), in the case of shear turbulence. Since 
the Navier-Stokes nonlinearity conserves both helicity and energy, one wonders whether there 
could exist situations characterized by an helicity cascade, analogous to the enstrophy cascade of 
two-dimensional turbulence || . 

All this neglects the presence of coherent structures and intermittency, which make an approach 
based on scaling and an hypothesis of homogeneous and isotropic inertial range questionable, even 
in the idealized case of spatially homogeneous large scales O, |lj, . The importance of hairpin 
vortices in wall turbulence [fl4[ [l5|| and even more, of plumes and effects at the boundary, in 
convective turbulence |ll| |17|, |18| is well known. Helicity, on the other hand, has long been 
suspected to play an important role in triggering intermittency in homogeneous turbulence |l^, . 

It should be mentioned, however, that with the exception of convective turbulence, intermit- 
tency and coherent structures seem to produce only minimal effects on energy spectra and transport 
coefficients. For this reason, they are not very interesting, when it comes to deriving turbulent 
models for engineering applications. 

The purpose of this paper is to study the behavior of the energy spectrum in the three examples 
of inhomogeneous turbulence listed. Dimensional analysis must necessarily be supplemented by 
information on the dynamics, to produce acceptable results. In wall turbulence, this will take the 
form of hypotheses on the distribution of vortices and on the way they are generated. In the case of 
convective turbulence, an analysis of the energy budget in the Navier-Stokes and the temperature 
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equations becomes necessary to verify the realizability of different scaling hypotheses. In helical 
turbulence, the same task is realized by means of closure analysis. 

In the next section, wall turbulence is analyzed, assuming that at any height, a Kolmogorov 
cascade is generated with integral scale equal to the height in exam. In section III, an analysis 
of the various possibilities for scaling in "homogeneous" convective turbulence is carried on, using 
also results from simulations of a GOY model of the type introduced by Jensen et Al. pi] , plus 
buoyancy couplings. The scaling predicted by Bolgiano |2^, ^3|, |], in particular, is taken under 
exam. Section IV is devoted to an analysis of helical turbulence using an EDQNM closure J24|, [To| . 
Section V contains the conclusions. 



II. Wall turbulence 

To fix the ideas imagine a turbulent flow parallel to a horizontal plane, characterized by a height 
5 and a stress at the surface (for unitary fluid density) v%. If the Reynolds number Re — ^ is 
very large, with v the fluid viscosity, the mean velocity V will obey to a very good degree of 
approximation the logarithmic profile law: jj| 

7x(x3) = ^log(^); r «X3«<5, (1) 

where k and (3 are dimensionless constants depending on the roughness of the wall, ro = v/v* is 
the flow inner length and X3 is the distance from the wall. The law of the wall, Eqn. (1), does not 
provide informations about turbulent fluctuations, beyond what could be obtained from a mixing 
length approximation; indicating by v the turbulent velocities: 

< V1V3 >= ~v T - — ; v T = lv T , (2) 
0x3 

where the eddy size I, the characteristic turbulent velocity vt and the eddy viscosity vt all depend 
on the height X3. In a scale invariant situation, using Eqns. (1) and (2): 

l(xa) ~ X3 and vt ~ «*, (3) 

which means simply that gradients of strength x 3 ~ 1 v* at scale X3 lead to vortices of size X3 and char- 
acteristic velocity u*. If one imagines that these vortices generate Kolmogorov cascades, spatially 
localized in height, some idea on the behavior of the energy spectrum can be obtained. 

This idea is not totally unreasonable, since, during the time it takes to the energy to be 
transferred to the viscous range, that is of the order of an integral time, eddies will move at most 
by a distance of the order of an integral length ~ X3 . 

Consider then the following picture (see Fig. 1): At a given height X3, eddies of size Iq ~ X3 
are generated by instability of the mean flow. These "mother eddies" split into "daughter eddies" 
of size lo a < Iq producing a Kolmogorov cascade; the index a indicates the point in the cascade 
and is the logarithm of the ratio of the size of the daughter eddy to that of the original mother 
eddy. To these however, there will be superimposed mother eddies generated above X3 and which 
will have therefore size l p > Iq ~ X3. Also these produce cascades superimposed with the original 
one, with daughter eddies of size l pa < l p . Indicate: 

l P * = e-% = e r~°X3 (4) 

It is clear that a > 0, but it is also true that p > 0. This last condition means that, at height X3, 
only mother eddies of size l p > X3 are present; smaller ones are generated at lower values of X3. 

In the presence of a Kolmogorov cascade, one will have that the typical ratio of the velocity 
inside daughter and mother eddies will be proportional to {l p<T /l p )^ exp(—r/l p(T ), with r ~ (353^) 31 
the viscous scale at height X3. One can then write for the velocity difference vi(x) = v(x+l) — v(x): 



v i( x ) ~ Y] (lp-) 3 exp(-r// Pi(Ti )ui(x,x^,p J ;,cr l ) (5) 



Pi 



3 



where p < R = log(S/x 3 ), and u(x, x',p, a) is the normalized velocity at position x, due to a 
vortex of type pa centered at x' (ui indicates finite difference with respect to the first argument). 




Figure f : Generation of Kolmogorov like cascade; for p > 0, l p is larger than the reference height 
x 3 at which the structure function S(l,x 3 ) is being measured. 

An assumption on the statistics of u(x, x', p, a) becomes therefore necessary. Spatial correla- 
tions in this model are assumed to arise from the finite extension of individual eddies, while distinct 
eddies are taken to be uncorrelated. This leads to the expression: 

< U (x, x'p, a)u(y, y', p' , a') >= v 2 6(p - p')S(a - a')6(x' - y')c(^f^) (6) 
with C(l~e\x — x'|) the normalized [C(0) = f] auto-correlation for a single vortex. In general, 

— (3— C ) 

vortices will be distributed with a density: n(x, p, a)l p „ "", with C, pa allowing for the inclusion 
of intermittency corrections in the model. However, if ro <C l pa -C l p -C 6, scale invariance implies 
n = n(lp 1 X3), with n(l) ~ f, while ( = for locally space filling (non-intermittent) turbulence. 
Using Eqn. (6), the two-point structure function can be computed explicitly: 

S(l,x 3 ) =v~ 2 <vf >~ dp daF(exp(\og(l/x 3 )-p + *))cxp(--a-^ye-^-rt) (7) 

where F(a) = 1 — C(a). Assuming smoothness of the velocity profile in an individual vortex and 
finiteness of its total energy leads to the conditions on F: 

F(a) = 0(a 2 ) (a < f) 

lim a 2 (F(a) - 1) = 0. (8) 

a— *oo 

For Re — ► oo, one can then obtain asymptotic expressions for S(l,x 3 ). For / <C x 3 <C S, the 
integrals are dominated by the contribution at p — and a = log(x 3 /l), i.e. the contribution from 
the eddies down the cascade generated at x 3 which have size I. In this way, Kolmogorov scaling 
arises: 

For I 3> x 3 the dominant contribution is from mother eddies with size x 3 <l p < min(7,<5). Due to 
the uniform distribution in p, this leads to a logarithmic form for the structure function: 

S(l,x 3 ) = -log(l/bx 3 ) + o((-y) x 3 ^l^S (9b) 
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S{1, x 3 ) = 77 log(6/bx 3 ) + e(x 3 /l) x 3 < 6 < I 



(9c) 



with e(x 3 /l) at most ©((a^//) 2 ). The coefficients a and & and the various higher order terms in 
Eqns. (9a-c), all depend on the detailed shape of the function F(a). The one-dimensional energy 

spectrum Q is: Ek(x 3 ) — J^(S(oo,x 3 ) — S(l, x 3 )) e lkl dl ~ dS ^ k Qk m tne inertial range. 

Hence, one has a standard k~i scaling for I <C x 3 and a fc _1 scaling when £3 <c I <C S. A definite 
expression for the energy spectrum can be obtained fixing the form of the auto-correlation. Taking 
C(a) = exp(— a 2 ), one obtains: 



E k {x 3 ) = n^v 2 x 3 



dp / da exp I p 



{kx 3 y 



a 2(p-<r) 



(10) 



A plot of this spectrum for different values of x 3 /S is shown in Fig. 2. A fit of wind tunnel data, 
taken from [SI, is shown in Fig. 3; the data correspond to values of the ratios: 2:3/15 ~ 0.01 and 
Tq/8 ~ 0.001, i.e. to the experiment with Reg = 7076 and y + — 28 illustrated in Fig. 1 of that 
reference. 

As discussed in vi samples in the range I ^> x 3 , velocities corresponding to vortices 

generated much above x 3 , while for I <C x 3l it samples the daughter eddies of size I generated in 
the cascade started at x 3 . In |27 , a fc" 1 scaling was obtained by means of a hierarchy of hairpin 
vortices of prescribed shape, but the transition to the k~i range was left out of the description. This 
phenomenon is explained here as a natural consequence of implementing a Kolmogorov cascade 
in an inhomogeneous turbulence setting; in all cases, hairpin vortices, and coherent structures in 
general, do not seem to be essential in obtaining such scaling behaviors. 
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Figure 2: One-dimensional energy spectra, for four different values of x 3 /S. a: x 3 /S = 0.001; b: 
X3/6 = 0.01; c: x 3 /5 = 0.1; d: x 3 /S = 0.5. In all cases: 6/r Q ~ 10 5 

The range / 3> x 3 is where turbulence ceases to be homogeneous and isotropic. Notice that 
the divergence of < v \ > as x 3 /S — > in Eqn. (9c) forces an implicit assumption of anisotropy 
in the model, in order to avoid inconsistency with Eqns (2). Thus, although < v\ > diverges as 
x 3 /8 — > 0, < v\v 3 > remains equal to w 2 ; this result can be obtained assuming that the velocity of 
vortices generated much above x 3 is almost parallel to the wall at £3. Hence, while < vf > receives 
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Figure 3: Fit of wind tunnel data for X3/8 ~ 0.01 and Vq/5 ~ 0.001. The curve is shifted to the 
left to overlap with the data, corresponding to a shift by ~ log(3), of the limits of integration in p, 
in Eqns. (7) and (10). Physically, this would correspond to a transition to anisotropy at I ~ 3x3. 

contributions from all vortices, up to size 5, < V1V3 > receives contribution only from vortices up 
to size £3. 

Although essentially kinematic, this model is dynamically consistent. In spite of the fact that all 
cascades overlap in space, the dominant interactions appear to be those among vortices of similar 
size, as in standard Kolmogorov theory, and belonging to the same cascade. A rough argument 
could be the following. The relevant strain over an eddy of size I, from a cascade starting at height 
X3, is produced by eddies of size V > I, This strain will be of the order of: vj' 1 (l'/x' 3 ) 1 ^ 3 with 
x' 3 the height of generation of the second cascade. If x' 3 < X3, however, the volume in which the 
interaction takes place will be reduced by x' 3 /x3 and the effective strain will be of the order of 
min(l, x' 3 /x3)v „V (I'/x^) 1 ^ 3 . The effective strain will be maximum then for X3 = x' 3 and I = I' 
corresponding to interaction with vortices of the same size and belonging to the same cascade. 

III. The case of "infinite space" convective turbulence 

Stratification in a fluid produces buoyant forces which couple velocity and temperature in the 
Navier-Stokes equation. In the limit of weak stratification in a large volume, one can adopt the 
Boussinesq approximation ||, in which the buoyant force in the Navier-Stokes equation, and the 
production term in the temperature equation, are linearized respectively in the temperature and 
the vertical velocity fluctuation. For a unitary density medium: 

(J^+v-v)v = VP + i/V 2 v-5e -1 e 3 (Ha) 

(^+v-v)0 = <tV 2 + 6V (116) 

Here, a is the molecular diffusivity, g is the gravitational acceleration, 9 and are the fluctuating 
and mean potential temperature, while 0' = The equations are linearly stable for 0' > 0; in 
this case, external forcing is necessary to achieve a stationary turbulent state. 
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In homogeneous isotropic turbulence, one derives Kolmogorov scaling [Q, using dimensional 
analysis with the quantities that are available: the scale I, the velocity difference vi and the mean 
(kinetic) energy dissipation e; this leads to the well known result: < vf >~ (eO 3 - l n the case of 
a stably stratified medium, in the presence of mechanical forcing, Bolgiano |2^, hypothesized 
an alternative situation, in which, potential energy transfer due to buoyancy forces, rather than 
kinetic energy transfer, governs turbulence dynamics. Dimensional reasoning with the quantities 
available in Eqns. (lla-b), leads then to the scalings ||: 

4 1 _ 2 

<v f >=Cvv (iy e fli; <vi6 l >=C vB (£)\l^; <0f>=CW(§) %J/*, (12) 
where eg is the dissipation of temperature fluctuations: 

€ 6 = <j\ve\ 2 . (13) 

It is interesting to carry on this dimensional reasoning, directly inside Eqns. (lla-b). After 
introducing the eddy turn-over frequency uii ~ I < vf > 5 and considering scales I much larger 
than the dissipation lengths for v and 6, we have from Eqn. (11a): 

uji < vf >~ ^ < v t 9i > and uji < vtfi ^ < Of >, (14) 

while, from Eqn. (lib): 

< viOi > ~ 6' < v\ ? > and lj ; < 9f >~ e e - 9' < > . (15) 

The question marks in Eqn. (15) indicate places in which the relation between terms is ambiguous. 
It appears that the ambiguity lies in the source term Q'v3 in Eqn. (lib). One sees immediately 
that the scaling described in Eqn. (12) is obtained when the term in 0' is negligible in Eqns. (lib) 
and (15). Introducing the Obukhov length: 

L„ = c Je'-i(|) 4 , (16) 

one realizes from Eqns. (12) and (15) that the condition of negligible 0' term is equivalent to 
I<L,. Thus, one has a source of temperature fluctuations at I < L„, whose energy is transferred 
to smaller scales by action of the convective term v • V0. These fluctuations provide the forcing 
for the velocity, through the buoyant term ^ff-, in Eqn. (11a). 

In order for a cascade of this form be present, the transfer of kinetic energy must be negligible, 
or equivalently, the frequency u>i (which is the strain felt by eddies at scale /) must dominate the 
one that would be produced in a Kolmogorov cascade; indicating with e the rate of kinetic energy 
production at the largest scales: 

uji ~ VL^L^^l^i > eH~i - v Lt L~*l~i (17) 

which implies 1 > L t . Thus, the two conditions of negligible 0' in Eqn. (15) and dominant 
buoyant force in Eqn. (14) restrict the possibility of Bolgiano scaling at most to a finite range 
around L*. 

An alternative situation, in which kinetic energy flows towards large scales, has been considered 
in p8| . In this case, the condition provided by Eqn. (17) is not necessary anymore; however, the 
condition of negligible 0' still forces Bolgiano scaling, to the range I < L*. Now, this is the range 
in which, the equation for the velocity decouples from that for the temperature, and it is difficult 
to see a mechanism whereby buoyancy could modify the nonlinear interaction in such a way to 
invert the direction of the energy transfer. If one restricts to this range, and maintains a situations 
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of forward energy transfer, the decoupling forces the temperature to be advected like a passive 
scalar, which results in the well known Kolmogorov-Corrsin scaling [ p9| : 

< vf >= C vv (d)i; < vi6 t >= C v eeg eUi < Of >= Ceeeee~hi, (18) 

The only way to obtain a power law fluctuation spectra for I > remains then, that: 

< v t di >~ g- l @e ~ O'^eg (19) 

This would lead again to a spectrum for < vf > and < 8f >; in this case, however, there would 
be a privileged scale, the Obukhov length, which would fix the amplitude of the cross-correlation 
< vi$i >: 



< viOi > (L 

,,2 ^ „- a2 



(< «? >< >) Z 

This leads to the situation, of correlations between vi and Oi being the strongest at I ~ and 
decaying at larger scales. 

From these observations, it is clear that scaling behaviors should not be expected for I > L*. An 
interesting question is then how, a toy system simulating a large range of scales like a GOY model, 
would behave in the presence of buoyancy. GOY models (see and references therein) present 
a cascade of energy along a linear chain of coupled ordinary differential equations for the complex 
variables u n , which are the analog of the velocity of eddies at scales l n = 2~ n in real turbulence. 
These models have attracted great attention, due to the coincidence of the intermittent properties 
of the moments < \u n \ n > and those of the structure functions < vf > in real turbulence. Jensen 
et Al. plfl have derived a generalization to the case of a passive scalar advected by a turbulent 
velocity held. It is easy to include in their model the effect of buoyancy; the resulting set of 
equations reads (k n = l~ ): 

(d t - vk 2 n )u n = ik n (au n+ iu n+2 + bu n -iu n+ i + cu„_iu„_ 2 )* - aT n + /„ 

(d t - crkl)T n = ik n [e n (u n -iT n+ i + u n+1 T n _x) 
+g(u„-2T n -i + it„-iT n _ 2 ) + h{u n+1 T n+2 + u n+2 T n+ i)}* + /3u n (21) 
for n = 1, 2, ...N, where the parameters a, b, c, e, <?, h are given by: 

= 1; b = c = g = -e=-h = -~; (22) 

and variables u m and T m in the nonlinearities are set identically equal to zero for m < 1 and 
to > N. The choice here is opposite to that of (28), in the sense that, it has been preferred not to 
tamper with the nonlinearities, and to leave them in the same form as without buoyancy. However, 
also in this case, some interesting results are obtained. 

Equations (21-22) have been integrated numerically for both stable and unstable conditions. 
In the stable case, a constant forcing /„ = (1 + i) * 10 _3 5„4 has been used to sustain fluctuations. 

In the unstable case no external forcing was present, and the fluctuations organized in such 

a way to push the Obukhov scale — log 2 (e T 2 /3 4 a 4 ) towards the lowest available shell; a 
Kolmogorov cascade ensued in all cases (see Fig. 4). 

In the stable case, the presence of two additional parameters with which to play: the forcing 
amplitude and wave number, allowed better control of N*. When TV* was sufficiently large, neither 
Bolgiano scaling, nor the situation depicted in Eqn. (18) took place, rather, a combination of the 
two, with steep ~ overlapping T- and u-spectra, and an almost constant cross-correlation 
< u n T n > (see Fig. 5). 

An attempt has been carried also to generate a backward kinetic energy cascade, using small 
scales forcing (the forced shell was n = 20). A physical interpretation of such a condition could be 
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Figure 4: GOY model simulation of convective turbulence under unstable conditions; values of the 
parameters: a — 0.01 [3 = —0.2; v = a = 10~ 7 ; no external forcing, a: < \T n \ 2 >; b: < \u n \ 2 >; c: 
< |u„T n | >. The buoyancy terms become of the same order of the others only for k n < 1; 




Figure 5: GOY model simulation of convective turbulence under stable conditions; values of the 
parameters: a = [3 = 10.; v = a = 10~ 7 , /„ = 10~ 3 (1 + i)S n 4. a: < \u n T n \ >; b: < \u n \ 2 >; c: 
< \T n \ 2 >. The buoyancy terms are of the same order of the others, over the whole range of k n , 
down to the dissipation scale. 



9 



the presence of plumes generated elsewhere in the fluid, in some steep thermal layer. Using Eqns. 
(21-22), the outcome was, however, a fc~ 5 / 3 spectrum in all situations. 

One can compare the behavior of in the two stability situations, looking at the ratios of the 
fluctuation amplitudes of the buoyancy terms aT n and (3u n , to those of dtu n and dtT n . One sees 
in Fig. 6 how in the stable case, buoyancy remains dominant over the whole inertial range, while 
the nonlinearity dominates the dynamics in the unstable case. 




10" 1 10° 10 1 10 2 10 3 10 4 



Figure 6: Plots of the ratios: r\ = ^rf~^f T25 and r 2 — ^\q^t\^> vs ' ^n- a: ri > stable, b: r 2 ; stable, 
c: n; unstable, d: r 2 ; unstable. 

In nature, of course, things go differently. First of all, the largest available scale lo corresponds 
to the size of the system. In the unstable case, steep thermal boundary layers develop rapidly and 
an approximation of constant temperature gradient ceases to be applicable. Thus, the left portion 
of the spectra in Fig. 4 is not particularly meaningful, and the prediction that — > Iq, under 
unstable conditions should not be trusted. In fact, in convective atmospheric turbulence, one has 
L* < Iq, and a large scale, buoyancy dominated range, is indeed present ]3l|| . 

Same amount of difficulties occur in the treatment of stable environments. In this case the 
problem is the idea of a purely large scale mechanical forcing. It is well known that, in these 
conditions, the large scale modes in turbulent shear layers are stabilized by buoyancy; thus, the 
peak in the forcing is moved to I ~ L*, and a range like that of Fig. 4 is not generated [ pj| . 
However, a purely large scale forcing can be generated, if low frequency waves are present in the 
flows; in this case, a turbulent spectrum extending to the buoyancy dominated region I < 
becomes possible again (see e.g. p2[). 



IV. Possibility of a purely helical cascade 

The last situation that is taken into consideration is that of non-reflection invariant turbulence. 
It is well known that, beyond energy, the nonlinearity of the Navier-Stokes equation has a second 
global invariant, the total helicity: 

H = \J d 3 xv(x) • [V x v(x)] (23) 

In many ways, helicity is the counterpart in three dimensions of two-dimensional vorticity, and a 
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natural question to ask is whether three-dimensional turbulence may exhibit multiple cascades, as 
it happens in two dimensions. The ability of helicity to hinder energy transfer p3| , |l0| , in particular, 
suggests the possibility of an helicity cascade with no energy transfer, given appropriate conditions 
on the forcing. 

Helicity, however, has the peculiarity of being a non-positive defined pseudo-scalar. Lack of 
positive defmiteness implies, in particular, that any triad of interacting modes, can exchange 
helicity in an arbitrary way, thus providing a source of helicity transfer fluctuations. 

This effect turns out to be important in GOY models; even in the case of maximum injection of 
helicity, for a given energy injection rate, it can be shown that the amount of the GOY equivalent 
of helicity: H = y]„ (~k n ) n \u n \ 2 |P0| , which is produced by fluctuations, is much greater than the 
one coming from forcing |34|j . Since the variable u n in GOY models mimics in a surprising way 
the velocity inside individual scale l n eddies, this may be a serious indication on the impossibility 
of an helicity cascade. 

Anyway, such a cascade seems impossible, also in a purely "mean field" description, with an 
helicity transfer to small scales, assumed constant over the whole space. 

The standard sequence of arguments, leading to Kolmogorov scaling, can be carried on, assum- 
ing a constant helicity flux eh to small scales; indicating with Hi ~ l~ 1 vf, the content of helicity 
at scale I, one can then write: 

e H - luiHi ~ I vf = const. =>■ vi ~ e 3 H l3 (24) 
implying expressions for the energy and helicity spectra and for the eddy turn-over frequency: 

E k = Cl e 3 H k 3; H k = c 2 e 3 H k 3; uj k = c 3 e 3 H k* (25) 

It is possible to obtain energy and helicity balance equations using statistical closure, starting from 
the expression for the velocity correlation —< |«k w -kl >: 

ki hi 

2irUl 3 = k- 2 P l J{k)E k + k- 4 e ijl k l H k ; P^'(k) = 6 ij - —. (26) 

k z 

Lesieur has derived such balance equations within the EDQNM closure [^0). In this kind of closure 
p4j , [To), the third order correlations < vvv >, which enter the equation for U l k 3 , are approximated 
by: < vvv >~< v w v m v m > + < v m v (1) v m > + < v^v'-^v 1 -^ >, with v {1> the first order 
pcrturbative solution to a modified Navier-Stokes equation, with the viscous term vk 2 replaced 
by the eddy turn-over frequency uj k . The v (0) are taken uncorrelated, so that the resulting 4- 
point correlations split into products of 2-point correlations. Furthermore, the approximation: 
< K0)w-k(°)l >- u k exp(-o; fe |t|) is adopted. 

The EDQNM equations for the energy and helicity balance read therefore JlO[ ] : 

(4 + 2uk 2 )E k - / dpdq9 kpq [—b kpq E q (k 2 E p - p 2 E k ) 
\dt ) J Ak pq 

l ^c kpq H q {k 2 H p -p 2 H k )\ (27a) 

P 9 



and 

d_ 

St 



2vk 2 \li k = [ dpdq9 kpq [—b kpq E q (k 2 H p -p 2 H k ) 

> Ja„ pq 



'A 

-—c kpq H q (k 2 E p -p 2 E k )}, (276) 



where: 

Okpq = (uj k + uj p + LUq)' 1 , b kpq = - (xy + z 3 ) and c kpq = z(l - y 2 ) (28) 
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with Afc the domain in which k, p and q can be the lengths of the sides of a triangle, and x, y and 
z the cosines of the angles opposite to these sides. 

The condition e# = const in Eqns. (24-25) and conservation of helicity triad by triad guarantee 
that the spectra and frequencies of Eqn. (25) provide a stationary solution for Eqn. (27b) for any 
value of the coefficients c,-. The energy balance, which is given by Eqn. (27a) fixes instead, at 
stationarity, the ratio of the two coefficients c± and C2- Numerical integration of that equation 
leads then to the result: 

— ~ 3.316 (29) 
ci 

However, from the definition of helicity, one has: 

fc -1 |.Hk| = 2vrfc < |[ik X v k ] • v_ k | >< E k = 2nk 2 < \v k \ 2 > (30) 

which implies |c2 < ci. Thus, Eqn. (29) cannot be satisfied, and an helicity cascade of the type 
described by Eqns. (24-25), does not seem to be possible. 



V. Conclusions 

The aim of this paper was to obtain some information, on the effect of large scale flow inhomo- 
geneities on the form of the energy spectra in turbulent fluids. Some idealized situations have been 
studied by means of simplified models and closure analysis. The point in common in the three 
inhomogeneous turbulence situations considered, is that simple dimensional reasoning, either gives 
wrong answers, or does not lead to any answer at all. Of course this was something to be expected, 
and in a certain sense, there is nothing deep in this result. However, the practical consequences 
are important. 

The analysis carried on here clearly shows that a fc _1 range is a universal feature of mechanical 
turbulent layers, which is independent of the presence of coherent structures. If one is interested 
in diffusion in wall turbulence situations, a k^ 1 range at scales kx% < 1 clearly makes a difference, 
with respect to a k~s spectrum, extending down to the inverse of the boundary layer thickness. 
Particles at distance I > X3 will separate horizontally, in almost a ballistic way: l(t) ~ t (with 
logarithmic corrections), while Richardson law: la(t) ~ t% will dominate in the vertical direction. 
The modifications that would be produced in dispersion models, for situations in which turbulence 
is predominantly mechanical, are clearly worth investigating. 

The interest for the existence of Bolgiano scaling is more academic, although some application to 
turbulent boundary layers in stable environments, in the presence of forcing by low frequency waves, 
is possible p2j . This scaling has attracted some interest a few years ago to explain observations 
carried on in liquid helium convection experiments p5| . This approach has been criticized later by 
several authors |l8], |3(|. The result of the analysis carried on here, suggests analogous difficulties 
for the existence of Bolgiano scaling in an idealized situation of convective turbulence in an infinite 
volume. The alternative however, which is characterized by velocity and temperature k~3 spectra, 
has difficulties itself due the presence of a privileged scale, the Obukhov length, dominating the 
dynamics, which weakens the very concept of an inertial range. GOY model simulations suggest 
indeed that neither scaling should be observed, rather, under stable conditions, a k~ 2 spectrum for 
both temperature and velocity, with correlation between the two, vanishing at large scales, should 
develop. 

The possibility of helical turbulence has sparked recently some attention in people interested in 
turbulence control f37f . A k~i energy spectrum, associated with an helicity cascade, would imply a 
decrease in the energy dissipation of the order of i?e _1 with respect to the standard k~i situation. 
(For equal total turbulent energy in the two cases: en ~ eL^ 1 =>• e' ~ Id^H with en and e' 
respectively the helicity and energy dissipation for afc~5 situation, e the energy dissipation for the 
corresponding k~% , and L and Id respectively the integral and viscous scales). The impossibility 
of an helicity cascade, suggested by the EDQNM calculation carried on here, means simply that 
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a k~3 could not be obtained modifying the large scale forcing, and that action at all scales (or 
alternatively at all frequency) would be necessary. Hence, turbulence control by forcing the cascade 
to become helicity dominated, could be possible only using a feedback system acting at inertial 
range frequencies. 
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